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Abstract 



' ^ ■ This paper gives a complete description of the solutions of the one 



dimensional Ginzburg-Landau equations which model superconductiv- 
fj ' ity phenomena in infinite slabs. We investigate this problem over the 

entire range of physically important parameters: a the size of the slab, 
K the Ginzburg-Landau parameter, and Hq, the exterior magnetic field. 
We do extensive numerical computations using the software AUTO, 
and determine the number, symmetry and stability of solutions for all 
^»0 i values of the parameters. In particular, our experiments reveal the 

^Jh ■ existence of two key-points in parameter space which play a central 

role in the formation of the complicated patterns by means of bifur- 
cation phenomena. Our global description also allows us to separate 
^^ I the various physically important regimes, to classify previous results 

in each regime according to the values of the parameters and to derive 
2 I new open problems. In addition, our investigation provides new insight 

into the problem of differentiating between the types of superconduc- 
, , tors in terms of the parameters. 
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1 Introduction 

This paper gives a complete description of the solutions of the one dimen- 
sional Ginzburg-Landau equations which model superconductivity phenom- 
ena in infinite slabs. One of our main goals here is to determine the number, 
symmetry and stability of solutions. We give a complete numerical investiga- 
tion of these issues over the entire range of physically important parameters: 
a the size of the slab, k the Ginzburg-Landau parameter, and /iq, the ex- 
terior magnetic field. Our experiments, which are summarized in Figure 
^, reveal the existence of two key-points in parameter space which play a 
central role in the way superconductivity is nucleated. Indeed, according to 
the values of a and k, superconductivity appears either in the volume of the 
sample if the slab is thin, or on superconducting sheaths near the boundary 
if the slab is large, and in the form of two-dimensionnal vortex patterns if 
the slab is of intermediate size. This is why all values of a and k are of 
interest. As we shall see, sheaths and vortices can be described in terms 
of asymmetric solutions of the one dimensional problem, and each of the 
key-points gives rise to diverse pattern formation by means of various bi- 
furcation phenomena. In addition, our results provide new insight into the 
problem of differentiating between the types of superconductors in terms 
of the parameters. Finally, our global description allows us to understand 
better many interesting physical phenomena, to classify results previously 
obtained for each regime of the parameters and also to derive new open 
problems. In order to properly describe our results we begin with a brief 
summary of the model. 

The superconductivity of certain metals is characterized at very low tem- 
peratures by the loss of electrical resistance and the expulsion of the exterior 
magnetic field hg. In the model derived by Ginzburg and Landau in 1950 
(see [^] ) , the electromagnetic properties of the material are completely de- 
scribed by the magnetic potential vector A (h = curl A being the magnetic 
field) and the complex- valued order parameter ip. In fact, ^ is an averaged 
wave function of the superconducting electrons and its modulus corresponds 
to the density of superconducting carriers. When the sample is wholly nor- 
mal, 1^1 = and the magnetic field inside the material h is equal to the 
exterior magnetic field hg. On the other hand, when the sample is perfectly 
superconducting, I';/;! = 1 and the magnetic field h is identically 0. Further- 
more, in the Ginzburg-Landau theory, the state of the sample is completely 
determined by the minimum of an energy depending on ip and A. For a 



more precise description of the general theory, one may refer to [19|, |24], 



§, H, |§,|^]orto|T8|,||]. 
In the special case when the sample is an infinite slab of constant thick- 
ness, between the planes x = —a and x = a, it is usual to assume that both 
■0 and A are uniform in the y and z directions, and that the exterior mag- 
netic field is tangential to the slab, that is hg =(0,0, /iq)- A suitable gauge 
can then be chosen so that -0 = f{x) is a real function, and A = q{x)ey, 
where e^ is the unit vector along the y direction (see [^ for more details). 
In this case, the nondimensionalized form of the Ginzburg-Landau energy is 
given by: 

EM,Q) = l[{^f" + fV + lf'-f + iq'-hof) dx. (1.1) 

The nondimensionalized parameter k is called the Ginzburg-Landau param- 
eter. It is the ratio of A, the penetration depth of the magnetic field, to 
^, the coherence length, which is the characteristic length of variation of 
/. The value of k determines the type of superconductor according to the 
type of phase transition which takes place between the normal phase and 
the superconducting phase: k small describes what is known as a type I 
superconductor and k, large as a type II. More precisely, for a type I su- 
perconductor, there is a critical magnetic field he such that if /iq < /lo the 
material is entirely superconducting and the magnetic field is expelled from 
the sample apart from a boundary layer of size A. This is called the Meissner 
effect. If ho > he, superconductivity is destroyed and the material is in the 
normal state, that is / = and q' = ho. For a type II superconductor, the 
phase transition is different and there are two critical fields hc^ and /ic2 • ^or 
ho < hci , the exterior magnetic field is expelled from the sample and there 
is a Meissner effect as for type I superconductors. But as ho is increased 
above hc^, superconductivity is not destroyed straight away, since the su- 
perconducting and the normal phase coexist under the form of filaments or 
vortices: the vortex is a zone of diameter ^, at the center of which the order 
parameter / vanishes. As ho increases further, the vortices become more 
numerous until the critical value /icj is reached at which superconductivity 
is destroyed. For ho > hc2, there is no superconductivity and the material 
is in the normal state. The way superconductivity is nucleated is highly de- 
pendent on a and k, as we will see later in Section 4, where we will introduce 
a third critical field hc^ corresponding to nucleation of surface superconduc- 
tivity. We refer to Tinkham ||3^ for a detailed explanation. The vortex 
phenomena in superconductivity have been widely studied in the literature. 
See for instance P] or p3] and the references therein. The critical value of 
K usually given to separate type I and type II superconductors is k = l/\/2- 



In the next Sections, we will describe how this value is computed in the lim- 
iting case a = oo. However, as a is decreased from infinity, we will see that 
the demarcation between type I and II behaviours is no longer the constant 
K = l/v2. Instead, we find that there is a well defined curve in the (k, a) 
plane which separates the two types of behaviours. This result will be fully 
discussed in Section 3. Similarly in Section 4, we find an extra curve which 
separates different type I superconductors according to surface effects. 

For a mathematical analysis of the problem, it is natural to assume that 
/ G H^(—a, a) and q G H^(—a, a). It then follows from standard variational 
arguments that there exists a minimizer of E^^, and that the minimizer is a 
solution of 

^r = f{f + q'-l) in {-a,a), 

q" = qf in {-a, a), ^'^'^> 

q'{±a) = ho. 

Notice that / = and q{x) = hQ{x + e) is always a solution for any real e. 
From now on, we will call this a normal solution. Regularity properties of 
minimizers yield that either / is a normal solution, or / does not change 
sign, hence we will study the case / > 0. An easy calculation shows that 
the energy E,^ is zero along the normal solution. Thus a global minimizer 
cannot have positive energy. 

The aim of this paper is to give a complete description (number, sym- 
metry and stability) of the solutions of the system (GL) for which / > on 
[—a, a], according to the values of the parameters a, k and /iq. This work is 
based on numerical simulations obtained with AUTO, a software developed 
by Doedel et al. (see ||2^, |2^, [^]) which computes bifurcation diagrams 
for systems of ODE's. 

In order to properly describe our results, it is necessary that we first give 
a brief summary of previous studies of the (GL) system. This is done in the 
next section. Then, in Sections 3 and 4 we give the details of our numerical 
investigation for symmetric solutions (Section 3) and asymmetric solutions 
(Section 4), together with their physical implications and suggestions for 
future research. Our results are all summarized in Figure |9[ Section 5 is 
devoted to the stability analysis. 



2 Mathematical background 

Let us first recall the basic properties of solutions. 

Proposition 2.1 // {f,q) is a solution of (GL) and if f is not identically 
zero then 

(i) I/I <l in {-a, a), 

(a) q has a unique zero oq in {—a, a), q is increasing on {—a, a), q' is 
decreasing on {—a,ao) and increasing on (ao,a). 

(Hi) There exist xi and X2 with —a < xi < oq < X2 < a and xq G [a:;i,2;2] 
such that f is increasing on [—a, xi]U[x2, a] and decreasing on [xi, X2\, 
f is increasing on [— a, xq] and decreasing on [x{),a\. 

The proof of {i) and {ii) can be found for instance in [^] and of {Hi) in [||]. 

There are two types of physically important solutions of {GL): symmet- 
ric solutions and asymmetric solutions. We define a symmetric solution to 
be a solution of {GL) such that f > 0, f is even and q is odd on [—a, a]. 
Thus, a symmetric solution satisfies the following problem: 

1 f// _ ff f2 I ^2 
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f" = f{f + q'-l) in (0,a), 



/(0) = /3, /'(0)=0, 
q" = qf^ in (0, a 
q{0) = 0, q'{0) = a 



2 ;„ (n „\ {GLsym) 



for j3 G (0, 1) and a > 0. We need to choose a and /? such that f'{a) = 0. 
Then {f,q) will be a solution of {GL) with /iq = q'{(i)- Notice that /? is the 
amplitude of / for a symmetric solution. 

We define an asymmetric solution to be a solution of {GL) which satisfies 
/ > on [—a, a], yet which is not symmetric. That is /'(O) / or q{0) ^ 0. 

For the existence of symmetric solutions, Kwong ||3^ has proved the 
following important result. 



Theorem 2.2 (Kwong jgZj/j For each (i in (0, 1), there exists a unique a > 
such that the solution {f,q) of {GLsym) satisfies f'{a) = 0. Moreover, a 
is a continuous, decreasing function of (3, 

lim a{(3) > and lim a{f3) = 0. (2.2) 



For this choice of a{P), and the corresponding solution {f,q) of (GLsym), 
let h{(3) = q'{a). Then h is well-defined, continuous and 

lim h{l5) = hs>0 and lim h(f3) = 0. (2.3) 

(3^0 13^1 



In the following the notation kg of (2^) plays an important role in the 
discussion of our results. 

There are three possible behaviours of the curve /i(/3) defined in Kwong's 
Theorem, and these are shown in Figures ||, ^ and |^. Notice that instead of 
graphing h{f5) vs /?, we have put /3 on the vertical axis and h on the horizontal 
axis, so that we keep the convention originally adopted by Ginzburg [p^ . 
Figures 3 and 5]. 

• Figure ||: for /5 G (0, 1), h{(3) is a decreasing function of (3: if < /iq < 
hs, there is exactly one symmetric solution of {GL), and if Hq > kg, 
no such solution. 



• 



Figure |2|: for [5 G (0, 1), h{(3) is increasing until h reaches a maximum 
value h (here /i ~ 1.1), and then decreasing to as /3 goes to 1. Thus, 
for Hq < hs, (GL) has a unique symmetric solution, for kg < Hq < h, 
two symmetric solutions, and for ho > h no such solution. 

• Figure g: for /5 G (0,1), /i(/3) is decreasing until h reaches a local 
minimum value h (here h ~ 0.83), increasing until h reaches a local 
maximum value h (here h ~ 0.95), and then decreasing to as /3 
goes to 1. There h < hg < h. Thus, for /iq < h, (GL) has a unique 
symmetric solution, for h < ho < hg, three symmetric solutions, for 
hg < ho < h, two symmetric solutions, and for ho > h no symmetric 
solution. 

If /i < hg, then for /iq < h, (GL) has a unique symmetric solution, 
for h < ho < h, three symmetric solutions, for h < ho < hg, a unique 
symmetric solution, and for ho > hg no symmetric solution. 

Note that in Figure ^ and ^, there are points on the curve where the resultant 
slope is vertical, when h = h or h. We shall refer to such points as folds 
throughout the paper. 

These figures indicate (and it has been verified experimentally) that if 
the applied field ho is sufficiently strong, superconductivity is destroyed. 
This has been proved by Kwong 



Proposition 2.3 (Kwong [Sq]) For any a and k, there exists he such that 
for ho > he, the only solution of (GL) is the normal solution f = 0, q' = ho- 
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Figure 1: Curve /i(/3) for a = 0.5 and k, = 0.4 
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Figure 2: Curve /i(/3) for a = 3 and k = 0.3 
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Figure 3: Curve /i(/3) for a = 3 and k = 0.9 



A weaker result had been obtained earlier in |38] and p9| . 

As each of the Figures indicates, as one decreases ho from infinity, the 
material remains in the normal state until a critical value of ho is reached 
at which there is a bifurcation of nontrivial solutions from the normal state. 
So for small e > 0, a nontrivial curve (/(.,e), g(.,e), /i(e)) of solutions of 
(GL) starting from a normal solution {0,ho{x + e),/io) is sought, with the 
following asymptotic development: 

f{x,e) =efo{x)+e^fi{x)+o{e^) in H^{-a,a), 

q{x,e) =qo{x)+e^qi{x) + o{e'^) in H'^(-a,a), (2.4) 

h{e) = ho + e^hi + e^/i2 + o(e^), 

where qoix) = hQ{x + e). It is important to note that when e = 0, the branch 
gives rise to symmetric solutions and when e 7^ to asymmetric solutions. 

Furthermore, it can be proved (see [0] or |^ for more details) that 
for small e > 0, the energy of the bifurcated superconducting solution, 
Ei^{f{.,e),q{.,e)), has the same sign as hi. This leads us to the following 
definition. 



Definition 2.4 A curve of the form ( \2.4 ) is said to result from a supercrit- 
ical bifurcation if hi < 0. Then, for small e > 0, the bifurcated solutions 
have lower energy than the normal state. 

A bifurcation is said to be subcritical if hi > 0. Here the bifurcated 
solutions have larger energy than the normal state. 

This means that for Figures || and ^, the bifurcation is supercritical and for 
Figure |2| subcritical. 

We now summarize the main results previously obtained concerning 
these figures describing symmetric solutions. We note that Ginzburg [g^] 
had investigated the case k small and found that h{f3) behaves as in Figure || 
for small a and, as a is increased through a critical value, the graph of /i(/3) 
changes from Figure |l] to y. He explained the type of behaviour described by 
Figure]^ in terms of superheating and supercooling. More precisely, when ho 
is large, superconductivity does not occur and the material is in the normal 
state. As ho is decreased, the material stays in the normal state down to hg, 
even though there is a range of ho where the normal solution is only a local 
minimizer and the global minimizer is a superconducting solution (see [ p^ . 
Figure 2]). If /iq is decreased further, there is a jump in the maximum of / 
and the material becomes superconducting, the solution being given by the 
symmetric branch. In this case, hg is called the supercooling field. Now, on 
the contrary, start from ho = where the superconducting state (1,0) is the 



global minimizer and increase /iq. The material will remain superconduct- 
ing until h is reached, though for fields slightly less then h, it is only a local 
minimizer and the global minimizer is the normal solution. For ho above h, 
there is a jump in the maximum of / and the material reverts to the normal 
state. This is the superheating phenomenon. These two phenomena give 
rise to a hysteresis loop as described in |2^ or later in |2^. We will make 
further remarks about these stability properties in Section 5. 

Ginzburg and Landau |2^ had also noticed that if o large is fixed, then 
there occurs a symmetric supercritical bifurcation of superconducting solu- 
tions from the normal state as n is increased through a critical value, but 
had no special explanation for this, since at that time only superconductors 
with small k were known. However, as we shall see in the next Section, 
as K is increased through this critical value, our studies indicate that the 
behaviour seen in Figure |2| changes into that seen in Figure |3|. 



Chapman [14| has studied the case = 00, and showed that there is a 
change of bifurcation from subcritical to supercritical that takes place for 
K = l/\/2) which is the critical value between type I and II superconductors. 
Moreover, in [^], a linear stability analysis through the time dependent 
equations yield that the value k = l/\/2 is also the one for which stability 
of the normal solution switches. In |]l6[, he has investigated the behaviour 
of the superheating field h in the limit a = 00. 

Following these works which are mainly based on formal computations, 
Bolley and Helffer have extensively studied the phenomenon of bifurcation of 
solutions from the normal solution in 0, §], [§, |l|], Q, H, |l|. In par- 



ticular, they have given rigorous proofs of properties of bifurcating branches 
and asymptotic formula for the superheating and supercooling fields. Al- 
though their results are mainly local, they have made a first attempt in [^] 
to give a global stability picture of the solutions. We will mention the details 
of their contributions relevant to the values of the parameters in the course 
of the discussion, and indicate how their results contribute to the broader 
global picture. 

Recently, Hastings, Kwong and Troy ||2^ have further investigated the 
multiplicity of symmetric solutions. They have proved that if k < l/\/2j 
and a is sufficiently large, then the behaviour of h{(3) is described by Figure 
g, and if K > l/v2 and a sufficiently large, by Figure ^. 

Seydel |3^ [^] was the first to give numerical evidence that there is a 
range of parameters for which asymmetric solutions and multiple symmetric 
solutions coexist. 



One of the major goals of this paper is to extend the results described 
above and completely determine the multiplicity and stability of symmetric 
and asymmetric solutions for the entire range of values of k and a. In 
the next two sections, we present the detail of our numerical investigations 
for symmetric solutions (Section 3) and asymmetric solutions (Section 4), 
together with their physical implications and suggestions for future research. 
In particular in Section 4, we explain how asymmetric solutions of {GL) give 
rise to nucleation of sheaths and vortices. The last Section is devoted to 
stability properties. 

3 Symmetric solutions 

As we stated in the last section, one of the main goals of this paper is to gain 
a complete understanding of the important properties of symmetric solutions 
in terms of the parameters k and a. As we shall see, these global properties 
give us new insight into an important problem, that of demarcating the 
(a, k) plane into regions corresponding to type I and type II behaviours. 
The results of our numerical investigations concerning the existence and 
multiplicity of symmetric solutions are shown graphically in Figure ^. They 
indicate that the (o, k) plane is the union of three connected sets Si, S2 
and 53. In 5i, the behaviour of h{f3) of Figure |l] holds. Likewise, ^2 and ^3 
reflect the behaviour of h{f3) of Figures ^ and y respectively. In the following, 
we explain the results shown in Figure |^ and make several conjectures. 

Conjecture 3.1 There exist three continuous functions Ki{a), ^2(0) and 
^3(0) which separate the (a, k) plane into three connected regions Si, S2 
and 53, as shown in Figure 0. There exists exactly one point a* called 
the triple point (with approximate value 1.60) such that Ki(a*) = ^2(0*) = 
^3(0^*) = K* — 0.93. Moreover, Ki{a) is defined and monotone increasing on 
[\/5/2,a*], ^3(0) is defined and monotone increasing on [a*, 00) and ^2(0) 
is defined and monotone decreasing on [a=K,oo), 

lim ^1(0) = 0, lim ^2(0) = l/v2, lim K^ia) = +00. (3.1) 

• region 5i = {0 < a < \/5/2,k > 0} U {\/5/2 < a < a^:,K > 
Ki{a)} U {a > a*,K > ^3(0)}. If (a, n) is in Si, the behaviour of 
the corresponding curve h[(3) is given by Figure |J and there is at most 
one symmetric solution of {GL). 

• region S2 = {\/5/2 < a < a^,0 < k < Ki{a)} U {a > a=K,0 < k < 
K2{a)}. If {a, k) is in S2, the behaviour of the corresponding curve h{0) 
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is given by Figure [^ and there are at most two symmetric solutions of 
(GL). 

Let us call (3 the point where h(P) = h. Then as (a, k) tends to a point 
(a,Ki(a)) on the curve ni{a), {(3,h) tends to (0,/is). 

• region S3 = {a > a*,K2(o) < k, < ^3(0)}. // (a, k) is in Ss, the 
behaviour of the corresponding curve h{(3) is given by Figure |^ and 
there are at most three symmetric solutions of {GL). 

Let us call (5 the point where h{P) = h. Then as (a, k) tends to a point 
(0,^2(0)) on the curve ^2(0), {I3,h) tends to (0,/is). As {a,K) tends to 
a point {d,Ks{d)) on the curve ^3(0), {P,h) and {f3,h) have the same 
limit (P, h{(5)) with (3 £ (0, 1). Moreover, (3 = if and only if a = a^,. 

Note that in 6*2, hs is sometimes called the supercooling field and in ^2 
and 5*3, h refers to the superheating field. We will come back to this in 
Section 5, since this is linked to local stability. 
Our physical interpretation of the 3 regions is the following: 
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* Si corresponds to what is usually called thin films: when the size of the 
slab is sufficiently small, there is no difi'erence between type I and type II 
superconductors. The Meissner effect is observed, that is the magnetic field 
is expelled from the material and there is no superheating nor supercooling 
phenomenon. We refer to |33] for more details. In particular, our curves 



Ki{a) and ^3(0) give a quantitative estimate of how small a should be de- 
pending on K so that only one dimensionnal phenomena happen. 

* We claim that the curve K2{a) for a > a^ can be considered as the sepa- 
rating curve between type I and type II behaviours. Indeed, it is known (see 



Tinkham [^^ for instance for details) that when /iq is decreased from infinity, 
type I materials become superconducting through a first order phase tran- 
sition and type II materials through a second order phase transition. More 
precisely, type I materials supercool, that is remain normal until /iq = hg, 
where nucleation occurs, followed by a discontinuous and irreversible jump 
in the maximum of /. This is due to the fact that the transition between 
the normal phase and the superconducting phase is accompanied by positive 
energy (the superconducting solution has a positive energy, that is bigger 
than the energy of the normal solution) , which corresponds to a bifurcation 
curve of locally unstable solutions (subcritical bifurcation). On the con- 
trary, for type II materials, / varies continuously at the transition because 
the superconducting solution has a negative energy, that is lower than the 
energy of the normal solution, so that the bifurcated curve of superconduct- 
ing solutions is locally stable. Therefore, we can take as a definition of type 
I (resp. type II) behaviours the fact that the bifurcation from the normal 
solution is subcritical (resp. supercritical). Hence, the curve ^2(0) sepa- 
rates the two regimes for finite a > a^ and we have lima—^oo ^2(0) = l/\/2: 
region ^2 corresponds to type I behaviours, where the bifurcation is sub- 
critical, and region 6*3 to type II, where the bifurcation is supercritical. In 
the limiting case a = 00, the formal computations of Chapman [O] indicate 
that the bifurcation curve of symmetric solutions changes from subcritical 
to supercritical as k passes through l/v2 from below. 



Below, we give further discussions and interpretations of the regions and 
the curves depicted in Figure §. 

3.1 Ki{a) 

For K = 0, Ginzburg |2|] obtained the critical value oq = \/5/2: if a < oq, 
h{P) is decreasing (Figure |l|) so that no supercooling nor superheating is 
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possible. For a > uq, h{j3) has the behaviour described by Figure ^ and 
displays both supercooling and superheating phenomena. 

Our study indicates that as (a, k) passes through Ki(a), the bifurcation 
at hs changes from supercritical in Si to subcritical in 52. Along the curve 
Ki(a), the function h{f3) is decreasing and we believe that /i"(0) = 0, or 
equivalently, the term /ii = in the development of h in (|2.4| ) . This means 
that as (a, k) enters 5*2, the fold in the h{j3) curve is going to appear near 
/3 = 0. The fact that /i"(0) = on the curve Ki{a) is directly linked to 
Ginzburg's observation that at the transition, the specific heat (which can 
be related to the inverse of a derivative of the energy) discontinuity is infinite. 
Further insight into the nature of this bifurcation could be obtained if one 
were to analyze a formula indicating the direction of bifurcation. Such a 
computation is directly linked to the evaluation of the term hi in (|2.4| ) . We 
refer to Millman-Keller |^l| for a formula for hi. 

3.2 K2{a) 

This is the limiting curve separating type I and type II behaviours as we 
have explained earlier. Our study indicates that as (o, n) passes from ^2 
to ^3 across ^2(0), there is a change from subcritical bifurcation in S2 to 
supercritical in S-^. Moreover, along the curve ^2(0), for a> a^ the function 
h{(5) has the graph indicated in Figure |2[ Additionnally /i"(0) = 0, since 
h'{(3) changes sign near /? = 0. This is similar to what is happening on 
Ki{a). Again by computing hi exactly as described above for Ki(a), one can 
obtain further mathematical results about the nature of the bifurcation and 
the location of K2{a)- 



The following additional information has been obtained in |13| about the 
asymptotic behaviour of the curve when a is large: 

Theorem 3.2 (B alley- Heljfer ^^) There exists a constant kq^q) > l/\/2 
such that for k > Ko{a) and for ho the eigenvalue of a spectral problem, then 
the bifurcated solution l[2.^) bifurcating from the normal solution (0, h^x, /iq) 
satisfies h{e) < h^. Moreover, there exists a constant C such that K,o{a) can 
be chosen such that 

—j=z < Ko{a) = —j=z + 0{exp{—Cc?y) when a — > cx). 

3.3 K3(a) 

Our study indicates that as (a, k) passes across ^3(0), the two folds at (/?, h) 
and {f5,h) which exist for a solution in ^3 tend to the same limit {(3,h). 

13 



As a result, the function /i(/3) is decreasing along ^3(0), but there is an 
interior point j3 with infinite slope, that is h"{p) = 0: this is where the 
curve gains its S shape as (a, k) enters 5*3 (see Figure ^). At the left end of 
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Figure 5: Part of the curve h{P) for a = 2 and k = 1. 

^3(0), that is at the triple point (a^,, n^,), we have /3 = 0. It is an interesting 
open mathematical problem to define the curve ^3(0) by means of a spectral 
problem in terms of h"{(3) = for some /? G [0, 1). 



3.4 The triple point 

As described earlier, our numerical experiments indicate the existence of a 
unique point a^, ~ 1.60 for which Ki(a=K) = ^2(0*) = ^3(0*) = ^*- Thus, 
in Figure |^ we see that this 'triple point' (a*,K*) provides the common 
intersection of the three curves. To our knowledge, the existence of this 
important point has not previously been reported in the literature. At the 
point itself, we observe that /i^ = /i* ~ 1.2, /i(/3) decreases on (0, 1) and we 
believe that /i"(0) = /i(*'')(0) = 0, that is /ii = /i2 = in (^. However, if 
we let (a, n) vary from the triple point along rays leading into the interior 
of ^i, S2 and 53, we expect to see entirely different behaviours. Along a ray 
leading into Si from the triple point, the bifurcation diagram for symmetric 
solutions is the same as that shown in Figure ||. Similarly, along rays leading 
into S2 and 5*3, the bifurcation diagrams for symmetric solutions are the 
same as those shown in Figure |2| and ^, respectively. 

Clearly, a careful bifurcation analysis needs to be done to expose the 
behaviour of solutions in a neighbourhood of the point a = a^,, k = k^,, 
h = h^. 
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3.5 Region Si 

Recall that the behaviour of h{P) is described by Figure |l[ 

3.5.1 K small 



This regime was originally studied by Ginzburg [25| in the case k = 0. 
He uses an approximate model in which / is held constant, i.e. / = /?. 
With this approximation, he observes that for a < v5/2, the curve h{/3) is 
decreasing, as described in Figure |l|. He computes hg = yS/a and claims 
that the approximate model remains valid in describing the properties of 
the solutions for small n. 

For the approximate model studied by Ginzburg, his computations have 
been made rigorous by Bolley-Helffer |9|- |l0|] . 

For the full system (GL), there is only the following local result con- 
cerning the bifurcation of symmetric solutions from the normal solution. 



Theorem 3.3 (Bolley-Helffer pq]) For any rj > 0, there exists Cq such that 
for Ka < Co, \/5/2 — a>r] and Hq the eigenvalue of a spectral problem, then 



the bifurcated solution (2.4) starting from the normal solution (0,hox,ho) 



satisfies h{e) < /iq- In particular, the bifurcation is supercritical. 
Rigorous treatment of the global problem remains open. 

3.5.2 K large 

Recently, this regime has received an increasing amount of attention since 
high temperature superconductors typically have large k. For mathematical 
results, see Berestycki, Bonnet, Chapman [Q] or Aftalion [y. 

Conjecture 3.4 Let a be fixed. When k tends to cxd, the limiting behaviour 
of the curve h{f3) is the union of a horizontal line at (3 = \ from h = to 
hQ^a), and a curve joining (/iQ(a),l) to {k,0), where /io(a) is given by 



a = j\u{hf-^{l-u'f 



As explained in Q, the limiting problem when k tends to cxd is 



/2 = (l-g2)i in (0,a) 



q" = qil - q^n\,\<i in (0,a), {GL^) 

q'{±a) = ho. 
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It is proved in Q that there is a unique solution (/oc'Zoo) of (GLoo) for 
ho < /iq and \qoo\ < 1 so that |/oo| > 0. This solution is symmetric, i.e. 
q{0) = and /(O) = 1. The critical field h^ is determined by the extra 
condition q{a) = 1. Moreover, any minimizer of £"« converges to (/ooQoo) 
as K tends to oo. Thus, it is natural to expect that the minimizer for finite 
K is symmetric and gives the horizontal line of the h{P) diagram. 

Note that lima_,oo ^o('^) ~ l/v2, which will be the same limit later in 
53 for a = 00, K large. 

Further study is needed to resolve the behaviour of solutions in this 
regime. 

3.6 Region 5*2 

Recall that the behaviour of h{/3) is described by Figure |2|. 

3.6.1 K small 

This regime was also studied by Ginzburg |25| with the approximate model 



where / is constant. He observes that for a > V^/2, the curve /i(/3) has 
the same behaviour as the one described by Figure ^. Again, only the 
computations for the approximate model have been made rigorous by Bolley- 



Helffer |11]. In particular, they prove that, for small k, the curve /i(/3) has 
a local maximum h, for which they give asymptotic formula consistent with 
physicists 'results . 

The only result to date for the full system {GL) concerns the bifurcation 
of symmetric solutions from the normal solution: 

Theorem 3.5 (BoUey-Helffer Ji^j For any r] > 0, there exists Gq such that 
for KG < Go, a — V^/2 > r/ and ho the eigenvalue of a spectral problem, then 



the bifurcated solution (2.4) starting from the normal solution (0, /iqx, /iq) 



satisfies h{e) > ho. In particular, the bifurcation is subcritical. 

3.6.2 a large 

Here the first rigorous multiplicity result is the following: 

Theorem 3.6 (Hastings, Kwong & Troy, W\U Let k < l/\/2 be fixed. 
Then for a large enough, there is a j3i in (1 — l/(Ka), 1) such that h{j3i) > hg. 
In particular, this means that S'2 7^ and that there are values of ho such 
that {GL) has two symmetric solutions. 
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From |27|, we also have the fohowing hmit, which was originahy predicted 
by Ginzburg-Landau |2q]: 



Proposition 3.7 (Hastings, Kwong & Troy, f^^) \vca.a^oo hg = k. 
The only bifurcation result available at the moment is the following: 



Theorem 3.8 (BoUey-Helffer jjq/j For any r] > 0, there exists Cq such that 
for na > Cq, 1/v2 — k > r/ and /iq the eigenvalue of a spectral problem, then 



the bifurcated solution (2.4) starting from the normal solution {O,hox,ho) 
satisfies h{e) > /iq- In particular, the bifurcation is subcritical. 

For the limiting behaviour of the curve h{P) when a is large, our exper- 
iments lead us to the following. 



Conjecture 3.9 Let k < 1/v 2 be fixed. When a tends to cx), the limiting 
behaviour of the curve h{(5) is the union of a horizontal line at (3 = \ from 
h = to /loo > 1/^2; o-iT-d a curve joining {h = Xj \pl, /3 = 1) to {h = k, (3 = 
0). 

The schematic limiting shape of the diagram can be found in Figure ^ 
or 14 at the end of Section 4. The curve 'Sym' corresponds to symmetric 
solutions. It means that the fold is getting close to the line (3 = 1 and the 



knee is getting very thin, as illustrated by Figure 12 in Section 4, where 
the extra curve corresponding to asymmetric solutions has been added. To 
help resolve this conjecture, one should investigate two relevant boundary 
value problems which we now describe: the first one is for the curve from 
(l/\/2, 1) to (k,0) in Figure 14. It is derived by letting a tend to oo in 



1 f// 


= f{f + q' 


-1) 


/(O) 


= /?, 


/'(O) = 


= 0, 


q" = 


qP 


in (0, 


do). 


QiO) 


= 0, 


q'{0) = 


a. 



m 



(0,oo) 



(3.2) 



Conjecture 3.10 Let k < l/v2 be fixed. For each (3 in (0, 1), there exists 

a unique a such that the solution of i\3.^ ) satisfies lim^^oo f{x) = 0. We 

define h{(3) = lim^^oo <l'{x). Then h{(3) is a continuous increasing function 

of (3, 

1 
lim h(f3) = K and lim hid) = —;=. (3.3) 

^^0 ^ ' /3^1 v^ 
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Theorem 3.11 Let k = l/v2- For each f3 in (0, 1), there exists an a such 
that the solution of f |3.^ j satisfies liniaj^oo f{x) = 0. Moreover, 



lim q'ix) = 1/V2 V/3G (0,1). 



This is proved by Chapman in [|14[ using the special decouphng of (GL) for 
K = l/v2 mentionned in [^: 

/ S = f' + fq/V2, 

\ F = q' + {f-l)/V2. 



(3.4) 



Chapman observes in |14] that a symmetric solution of (GL) leads to (>S'(0), 
F{0)) = (0, 0). It follows from uniqueness of solutions of a first order system 
that {S, F) = (0, 0) and therefore the full (GL) system reduces to 

/ /' = -fq/V2, 
1 q' = {l-p)/V2. 



(3.5) 



(3.6) 



The second BVP relevant to Conjecture |3|£ corresponds to the horizontal 
line /3 = 1 in the bifurcation diagram of Figure IJ. It comes from a boundary 
layer analysis: the extremal part of the slab is sent to and the symmetry 
plane to cx), so that we have 

f ^f" = fif + q'-l) in (0,cx)), 
/'(O) = 0, lim.^oo fix) = 1, 
q" = qf^ in (0,cx)), 
q'{0) = h, lima;^oo 9(2;) = 0. 

This problem and the boundary layer analysis were originally derived by 
Ginzburg [^5[, and |£^, Figure 5] leads to the following conjecture: 

Conjecture 3.12 Let k > be fixed. There exists /iqo > l/\/2 such that 
(3.6) has a unique solution for h < l/-v/2 and two solutions for l/\/2 < h < 

This explains why the curve going down from /5 = 1 to /3 = in Figure ^ 
starts at h = l/v2- For this conjecture, partial results have been provided 
by Bolley-Helffer: 

Theorem 3.13 (Bolley-Helffer ^l^j) For any k > 0, there exists /ioo('^) 
such that there is no solution of (jS.dj) if h > hoo{K), and at least a solution 
if h < hoc . Moreover, 



V2 
4 ■ 
This formula was first derived by De Gennes. 



limK/i^(K) 



(3.7) 
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3.7 Region S3 

Recall that the behaviour of h{P) is described by Figure |3[ 

3.7.1 a large 

Here the only rigorous multiplicity result is the following: 

Theorem 3.14 (Hastings, Kwong & Troy, ^^) Let k > l/\/2 be fixed. 
Then for a large enough, there are /?i and (32 in (1 — l/(fi;a),l) such that 
(^2 > /3ij h{Pi) < kg and h{P2) > h{Pi). In particular, this means that the 
set 53 7^ and that there are values of Hq such that (GL) has at least 3 
solutions. 

Conjecture 3.15 Let k > l/\/2 be fixed. When a tends to cxd, the limiting 
behaviour of the curve h{(3) is the union of a horizontal line at P = 1 from 
h = to /loo > l/v2; o-nd a curve joining {h = l/v2, /3 = 1) to {h = k, P = 
0). 

Notice that the bifurcation curve for a large has the same limiting shape as 
in S2, but this time the curve from (l/\/2, 1) to (k, 0) has negative slope. 



This is illustrated by Figure 15 below with the curve 'Sym' corresponding 
to symmetric solutions. The part of the diagram corresponding to the line 
/3 = 1 is described by the same limiting problem as in 52. Thus, the field 



/loo is the value derived from Conjecture |3.12 . 



Conjecture 3.16 Let k > l/\/2 be fixed. For each /3 in (0, 1), there exists 
a unique a such that the solution of l\3.^ ) satisfies lima;^oo f'{x) = 0. We 
define h{P) = lim^^oo Q'ix). Then /i(/3) is a continuous decreasing function 
of P, and 

lim /i(/9) = K and lim /i(/3) = —^. (3-8) 

/3^0 /3^1 \/2 

3.7.2 K large, a = 00 

The asymptotic behaviour of the upper bound for the superheating field 
hooin) in 5*3 has been proved in the following: 



Theorem 3.17 (Bolley-Heljfer \1^]) The field /ioo(k) defined in Theorem 



3.13[ satisfies 

liin^/ioo(K) = ^. (3.9) 

V2 



K— >00 
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This proves part of the formal asymptotic development derived by Chapman 

[16|: hoc = l/\/2 + 0.33fi;-^/3 + o{k-^/^) as k tends to oo. Notice that 



this means that the 'superheating part' of the curve is getting small (recall 
Conjecture 3.16| ) and is consistent with the fact that we reach the boundary 



of 5i. We find the same limit as lima^oo ^o(^)- 

3.7.3 Ka large 

The only result so far concerning the bifurcation of symmetric solutions in 
this regime is the following 



Theorem 3.18 (BoUey-Helffer \19iJ) For any 77 > 0, there exists Cq such 
that for Ka > Co, k — l/v2 > rj and Hq the eigenvalue of a spectral prob- 
lem, then the bifurcated solution ( \2.^ starting from the normal solution 
(0,/iox,/io) satisfies h{e) < Hq. In particular, the bifurcation is supercritical. 

4 Asymmetric solutions 

Recall from Section 2 that (/, q) is an asymmetric solution of (GL) if it 
is not symmetric. Notice that if (/(x), q{x)) is a solution, then so is 
{f{—x),q{—x)), but in the following we will count this as one solution in- 
stead of two and assume that the maximum of / is reached in [—a, 0]. Below, 
we summarize the results of our numerical investigation of asymmetric so- 
lutions shown in Figures y and 0: the vertical axis has been labeled ||/||. 
For the symmetric branch, ||/|| refers to /(O) = (3 as in Figures ||, |2| and y. 
For the asymmetric branch, ||/|| = ||/||oo) which is attained at some point 
Xj3 different from zero. For simplicity, we will sometimes refer to the com- 
plete bifurcation diagram as /i(/3). In Figures ^ and |^, the curve /i(/3) for 
symmetric solutions is the typical curve in 82- On these curves, there is a 
point {Pb, /ife) denoted by a square: it is the branching point of the branch 
of asymmetric solutions. More precisely, for ||/||oo close to Pb, the asymmet- 
ric solution is nearly symmetric. The branch of asymmetric solutions leads 
from the branch of symmetric solutions to the branch of normal solutions. 
If the branch of asymmetric solutions crosses the branch of symmetric so- 
lutions without any square indicated, it just means that for this special ho, 
(GL) has a symmetric solution and an asymmetric solution having the same 
maximum value of /. We have observed the following: 

Conjecture 4.1 Let xp be the point where f reaches its maximum. Then 
Xfs is a decreasing function of (3 and lim^__^^j^ xp = 0. 
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Figure 6: Bifurcation diagrams for a = 3 and k = 0.35. 
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Figure 7: Bifurcation diagrams for a = 3 and n = 0.5. 
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We now describe the two possible behaviours of the asymmetric curve h{f3). 

• Figure ^ for /? G (0, /3fe), /i(/3) is an increasing function of (3 and 
limp^oh^P) = has, Iinii3^/3^h{f3) = /i^, with hs < has < h < h. If 
ho € {has,hh), there is a unique asymmetric solution of (GL) and if 
ho < has or ho > hi), no such solution. 

• Figure |^: for (3 £ (0, /?{,), /i(/3) is a decreasing function of /3 until /i 
reaches a minimum value h^g and then increasing to hb as /? goes to 
Pb- thus, for ho G (Z^as' ^fe)> there are two asymmetric solution of (GL), 
for ho G {hb,has), a unique asymmetric solution and for ho < has oi' 
ho > hb, no such solution. 

This implies in particular that for Figure ^ the bifurcation of asymmetric 
solutions from the normal solution is subcritical and for Figure ^ supercrit- 
ical. 

In Figure |8|, we have computed two new curves, Ki{a) and ^5(0) which 
play an essential role in the study of existence and multiplicity of asymmetric 
solutions. 

Conjecture 4.2 There exist two continuous functions ^4(0) and K^{a) sep- 
arating the (a, k) plane into three connected regions Ao, A\ and A^. There 
exists exactly one point a* called the quintuple point (with approximate value 
1.23) such that Ki{a*) = K5{a*) = l/\/2 and l/\/2 = Ki{a*). Moreover, 
Ki{a) = C/a with C approximately equal to 0.90, and K^{a) is defined and 
monotone decreasing on [a*, 00) with Ivoia^oo i^b{o) = f^as — 0.4. 

• region ^0 = {^^ > 0, < k < ^4(0)}. There are no asymmetric 
solutions of {GL). Thus in the h{(3) plane, we only see curves as 
shown in Figure |^ and |^. 

• region A\ = {a* < 0,^4(0) < n < n^{a)}. If (a, k) is in A\, the 
bifurcation diagram for asymmetric solutions is given by Figure & it 
has no fold and there is at most one asymmetric solution of (GL). 

Moreover, when (o, k) tends a point (a, ^4(0)) on the curve ^4(0), then 
has cind hb tend to hs. When (a, k) tends to some (a, ^5(0)), then has 
tends to hs. 

• region A2 = {0 < a < a* ,k > ^4(0)} U {a > a* ,k > ^5(0)}. // (a, k) 
is in A2, the bifurcation diagram for asymmetric solutions is given by 
Figure ^, it has a unique fold and there are at most two asymmetric 
solutions of{GL). 
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Moreover, when (a, k) tends to a point (a,K4(a)) on K4,{a), then has 
and hfj tend to hg. When (a, k) tends to a point (0,^5(0)) on K5(a), 
then has tends to hg and has ^^ ^b- 

Of particular physical interest is the following: 

Conjecture 4.3 Let a> 0. If n > K4^{a) then has > ^s- 

Thus, as ho decreases from infinity, the material is in the normal state 
until /iQ reaches the value has where the bifurcation of asymmetric solutions 
occurs. If a > a* and k. > ^5(0), or if a < a* and k, > K4(a), then the 
bifurcating asymmetric branch is supercritical and plays an important role 
in the onset of superconductivity in the material. In the following we explain 



the implications of Conjecture 4.3. In the physics literature, has is defined 
as /ic3 and hs as hc2- It is well known to physicists that when k is fixed 
bigger than 0.4 and a is taken large enough (which corresponds to k being 
above max(K4(a), ^5(0))), when the exterior magnetic field /iq is decreased 
from infinity, superconductivity is not nucleated first in the volume of the 
sample, which would give rise to symmetric solutions, but rather in a sheath 
near the surface, due to the existence of asymmetric solutions. This is called 
surface superconductivity. If the slab is very thick, the two surface solutions 
f{x) and /(— x) do not interact. In this region, superconductivity is first 
nucleated in surface layers of size 1/k near the boundaries, and the middle 
part of the material is normal. Now, if the slab is of intermediate size, the 
solutions /(x) and f{—x) interfere to create vortices. Indeed, recall that 
the original Ginzburg-Landau energy is gauge invariant so that a solution 
{f{x),q{x)) has the same energy as {e^'^'^y f{x),q{x) + c) for any constant 
c. Thus, when the sample is not too large, and h decreases below has, the 
linear combination of the asymmetric solutions f{x) and f{—x) create two 
dimensionnal vortices along the mid-plane x = 0. This is reflected in the 
formula: 

V^ = cos ky{f{x) + /(-x)) + i sin ky{f{x) - /(-x)) (4.1) 

Further discussion of vortex formation and the details of the derivation of 



formula (|4.1| ) are given in Tinkham |37|. In particular, the distance between 



two vortices along the x = plane is proportionnal to 1/x^. Indeed, when 
moving up the asymmetric branch /i(/3), we find that xp tends to zero as [3 
approaches the branching point on the symmetric branch. This means that 
the distance between the two vortices tends to infinity. At the limit, the 
material is perfectly superconducting because / is symmetric. 
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In fact, as described in |3^], hs and has can both be measured: has 
by nucleation of superconductivity as we have seen and hs by magnetic 
transition, which needs bulk superconductivity to take place. 

Boeck and Chapman Q have studied surface superconductivity in detail 
and they have determined the regimes of (a, n) where the asymmetric solu- 
tion gives rise to a surface sheath or to vortex solutions, through formula 

Now recall that the surface sheath is of size 1/k (that is of order of the 
coherence length ^), so it is consistent that a condition for surface super- 
conductivity to exist should be a > C/k. As mentionned by Tinkham p7| , 
p. 136], physicists had found that 2C ~ 1.81 but did not know how to predict 
the details of changeover of behaviour on this curve. We believe that our 
numerical results answer this open question. A curve similar to ^4(0) was 
also mentionned by Boeck and Chapman |^]. 



The global results of our numerical investigations are shown graphically 
in Figure ^ where for each region of interest, we refer to the relevant Fig- 
ure describing the bifurcation diagram of both symmetric and asymmetric 
solutions. 
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Figure 10: Bifurcation diagrams for a = 3 and k = 0.9. 
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Figure 11: Bifurcation diagrams for a = 1 and k = 1.2. 



Recall that a fold in the symmetric branch correponds to a superheat- 
ing phenomenon, a subcritical bifurcation at the highest nucleation field to 
a supercooling phenomenon and a supercritical bifurcation of asymmetric 
solution to surface superconductivity. In particular, when surface supercon- 
ductivity occurs there is no supercooling. Thus, our physical interpretation 
of our new regions in the light of the discussion given by Tinkham [ p7| is 
the following: 

* K < ^4(0) or K > Ki(a) and a < a* , that is Si D Aq (Figure 1): thin films 
with type I behaviour. There is no supercooling nor superheating. 

* K > K4(a) and k > Ki{a) and k > ^3(0), that is 5*1 n A2 (Figure 10): thin 
film with type II behaviour. There is no supercooling nor superheating, 
but there is surface superconductivity due to the presence of a supercritical 
branch of asymmetric solutions. 
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* K < K5(a) and k < Ki{a), that is S2 n Ai and S2 n Aq (Figures 2 and 6): 
type I behaviour. These superconductors supercool, that is the bifurcated 
branches from the normal solution are subcritical and superconductivity is 
nucleated in the bulk. 

* K5(a) < K < K2{a) and a> a*, that is 82^1 A2 (Figure 7): the material has 
type I characteristics because at hg the bifurcation of symmetric solutions is 
subcritical. However, the bifurcation of asymmetric solutions is supercriti- 
cal, and this results in the formation of surface sheaths of superconductivity. 



Thus, in this regime, Tinkham [37| has given the name Type I 1/2 to the 
superconductor . 

* K2{a) < K < Ks{a), that is S3 n A2 (Figure 9): type II behaviour, with 
superheating and surface superconductivity. 

Thus, we claim that the curve ^5(0) is the limiting curve separating type 
I behaviours and type I 1/2. 

Note that surface sheath effects can be inhibited if the superconducting 
surface is coated with a normal metal. In this situation, the boundary con- 
dition /'(iba) = needs to be changed and replaced by /'(±a) = —f{±a)/b 



as introduced by de Gennes |24]. This is why the type of superconductor is 



determined by the behaviour of the bifurcated branch at hg and not at has ■ 

To our knowledge, the global behaviour of the asymmetric bifurcation 
curves have never been mentionned in the literature. In light of the discus- 
sion given above the following open problems have physical importance. 

* prove how asymmetric solutions arise from bifurcation from symmetric 
solutions; 

* prove that as soon as the asymmetric curve arises from a subcritical bi- 
furcation from the normal solution, the curve has no fold; 

* prove that as soon as the asymmetric curve arises from a supercritical 
bifurcation from the normal solution, the curve has a unique fold. 

4.1 The quintuple point 

As mentionned earlier, our experiments reveal the existence of a second key- 
point, in addition to the triple point defined in Section 3. To our knowledge, 
the existence of this point has not previously been reported in the literature. 
For our discussion of the properties of this point, we refer the reader to 
Figure ^ First, we let a* be the unique value such that Ki{a*) = l/\/2. 
We have o* ~ 1.23. Then our investigation shows that (a*, l/v2) is on the 
curve K4(a) and is the end point of ^5(0). At the point itself, we observe 
that h{(3) is decreasing to hg = h* c^ 1.45 and there are no asymmetric 
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solutions. However, if we let (a, k) vary from this point, we expect to see 
5 different behaviours of the bifurcation diagram. Along a ray leading into 
SiHAq, the bifurcation diagram is the same as shown in Figure Q and there 
is no asymmetric solution. If (a, k) lies on a ray leading into Si Ci A2, the 
bifurcation diagram is the same as shown in Figure ^: the asymmetric 
solutions bifurcate from h* and the asymmetric branch has a fold, though 
it is not easy to see it on the picture. Next, let (a, k) lie on a ray leading 
into 52 n ^o; the bifurcation diagram is as in Figure Q so that the fold in 
the symmetric branch results from a subcritical bifurcation. If the ray leads 
into 52 n ^1, the bifurcation diagram is as in Figure |6| so that there is a 
fold in the symmetric branch resulting from a subcritical bifurcation at hs 
and there is also an asymmetric branch, with no fold, bifurcating from has- 
Both hs and has are very close to h* . Finally, if (a, k) lies in 52 H ^2, the 
diagram is as in Figure |^, this time the branch of asymmetric solutions has 
a fold. 



As we have seen in the discussion following Conjecture 4.3, the asymmet- 
ric solutions which appear play an important role in 2D vortex formation for 
intermediate values of a larger than a* . In fact, our investigation indicates 
that their existence is due to a complicated bifurcation phenomenon which 
occurs at the quintuple point. Thus, a complete mathematical investiga- 
tion of bifurcation phenomena occuring at the quintuple point needs to be 
carried out. 

4.2 K4{a) 

Along the curve ^4(0), there are only symmetric solutions. As (a, k) passes 
from ^0 to Ai or A2 across ^4(0), the asymmetric solutions bifurcate from 
f3 = 0, h = has and has tends to hs as (a, k) approaches a point on ^4(0). 

The nonexistence of asymmetric solutions bifurcating from the normal 
solution for na small was obtained by Bolley-Helffer [O]. 



Theorem 4.4 (Bolley-Helffer, pjjl) There is a constant Cq such that for 
KG < Co, there exists a unique nontrivial C°° curve {f{.,e),q{.,e),h{e)) 
of solutions of {GL) bifurcating from a normal solution (0,/io(a; + e),/io). 
These solutions bifurcate from the particular normal solutions {0,hQX,hQ) 
where /iq is the eigenvalue of a spectral problem, and therefore are symmetric 
solutions. 

The existence of asymmetric solutions bifurcating from the normal solu- 



tion for Ka large was obtained by Bolley-Helffer |13|. 



27 



Theorem 4.5 (Bolley-Heljjer, [T^) There exists a constant Ci such that 
for Ka > Ci, there exist exactly three nontrivial C°° curves {f {.,£), q{.,e), 
h{£)) of solutions of (GL) bifurcating from normal solutions (0, /io(x+e), /iq), 
where Hq is an eigenvalue of a spectral problem. One curve of solutions 
starts from the particular normal solutions (0, /iqx, Hq), and is a curve of 
symmetric solutions. The two others correspond to e = e(K, a) and e = 
—e{K, a) and are asymmetric solutions. 

These two Theorems only give local behaviour of the bifurcation curve near 
the normal solutions, but this scaling is consistent with the behaviour of 
Ki{a) which is C /a. 

A global treatment of the region below ^4(0) for the full Ginzburg- 
Landau system remains an important open problem. 

4.3 K^{a) 

As we have seen, K^{a) is an extremely important curve. Our study indi- 
cates that as (a, k) passes from Ai to A2 across ^5(0), there is a change of 
bifurcation for the asymmetric branch from subcritical in Ai to supercritical 
in A2 ■ This behaviour of the asymmetric branch is similar to what happens 
to the symmetric branch across Ki(a) and ^2(0). 

The following additionnal information has been obtained in [^] about 
the asymptotic behaviour of ^5(0) when a is large. 

Theorem 4.6 (Bolley-Helffer, pW) There exists a constant Sq approxi- 
mately equal to O.4 such that for any r] > there exists C2 such that for 
Ka > C2, \k — T,q\ > rj, e = e(«;, a) and ho the eigenvalue of a spectral prob- 
lem, then the solution of the form ( ^.4 ) bifurcating from the normal solution 
(0, ho{x + e), ho) satisfies (k — So)(/i(e) — /iq) < 0. In particular, if n < Sq, 
the bifurcation is subcritical, and if k > T,q the bifurcation is supercritical. 

Notice that Sq is what we call Kas- In lH^! ^0 is given in terms of two 



integrals. This value already appears in physics papers (see [32| for instance) 



and it is known to separate different behaviours of type I superconductors. 

4.4 a large 

The only rigorous result concerning the behaviour of the asymmetric branch 



is the following, which can be illustrated by Figure IC. 

Theorem 4.7 (Hastings and Troy J^) Let k > 1/V2.01 be fixed. Then 
for a large enough, there is a range of ho for which there is no symmetric 
solution and yet there is an asymmetric solution. 
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This Theorem imphes in particular that in the (a, k) range mentionned, we 
have has > h and has > hs- This gives a partial resolution to Conjecture 



4.3. The physical importance of this result will be discussed in Section 5. 



Figure |T^ illustrates the behaviour of the bifurcation diagram when a 
is large. This is the limiting case of the one shown in Figure y when a 
gets large. We see that the branching point of the asymmetric curve has 
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h 

Figure 12: Bifurcation diagrams for a = 10 and k = 0.3. 

the same limit as the fold of the symmetric curve. Thus, for the limiting 
behaviour of the asymmetric curve when a is large, our experiment lead us 
to the following. 

Conjecture 4.8 Let k he fixed. When a tends to oo, the limiting behaviour 
of the bifurcation curve of the asymmetric solutions is the union of a hor- 
izontal line at P = 1 from h = l/v2 to h = h^o (defined in Conjecture 
3.1^:) and a curve joining {h = l/\/2,(3 = 1) to {h = K/fi^,(3 = 0), where 



Hi is approximately 0. 59, and is defined as the minimum over a of the first 
eigenvalue of the Neumann problem of the harmonic oscillator in (— a,oo). 
In particular, 



lim hf, = /loo <ind lim hg 



1/V2. 



(4.2) 



This implies that k^s = //i/v2: when k < Kas, the curve from /3 = 
to /3 = 1 has positive slope resulting from a subcritical bifurcation (see 
Figure ^) and when k > Ka^, this curve has negative slope resulting from 
a supercritical bifurcation (see Figures 14 and |l5|). 

The value of /ij* and the corresponding spectral problem were derived by 
Bohey-Helffer: 
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Figure 13: Limiting bifurcation diagrams for a = co and n = 0.3. 
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Figure 14: Limiting bifurcation diagrams for a = oo and k = 0.5. 
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Figure 15: Limiting bifurcation diagrams for a = oc and k = 0.9. 
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Theorem 4.9 (Bolley-Heljjer, J^^) Let k be fixed. Then 

lim has = -K. (4.3) 

This is consistent with the experimental value found by St-James and De 
Gennes |2| which is 2.392k/ ^/2. 



For fixed k, Chapman, Howinson, McLeod and Ockendon p^] prove the 
existence and uniqueness of the solution of 



< 



' ^r' = f{f + q'-l) in (-00,00), 
lim^._^_oo fix) = 1, limj._^oo f{x) = 0, 
q" = qf in (-00,00), 



(4.4) 



lim^j^oo g(a:) = 0, lim^_^oo q''(a;) - ^ 



This gives the point P = 1, h = l/\/2 in the bifurcation diagram, but further 
study is needed to resolve Conjecture [4.8| . 

5 Stability of solutions 

We define (/, q) to be locally stable if it is a local minimizer of £"« in H^ and 
unstable otherwise. Local and global stability properties of both symmetric 
and asymmetric solutions have received recent attention in the literature. 
As a result of our investigations, we make the following conjecture, which is 
to be read in the light of Figure ^. 

Conjecture 5.1 In Si Ci Aq, the symmetric solution is a global minimizer 
when it exists, that is for < Hq < kg and the normal solution is a global 
minimizer for Hq > kg. 

In S2 n ^0 o,nd in S2 D Ai, the symmetric solution is locally stable for 
< ho < h, the normal solution is locally stable for h^ > hg, and in S2 D Ai 
the asymmetric solution is unstable. Moreover, there exists he £ {hg, h) such 
that the symmetric solution is a global minimizer for < ho < he and the 
normal solution is a global minimizer for ho > he. 

In Si n A2 the symmetric solution is locally stable for < /iq < hb, the 
asymmetric solution is locally stable for h^is < ^0 < ^as (^nd the normal 
solution is locally stable for /iq > hag. Moreover, there exists he G (has^hs) 
such that the symmetric solution is a global minimizer for < ho < he, the 
asymmetric solution is a global minimizer for he < ho < hag and the normal 
solution is a global minimizer for ho > hag ■ 
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In 5*2 n A2 the symmetric solution is locally stable for < Iiq < h, the 
asyvfimetric solution is locally stable for h^^ < Hq < has CLnd the normal 
solution is locally stable for /iq > has- Moreover, there exists he G {hs,h) 
such that the symmetric solution is a global minimizer for < hQ < he, the 
asymmetric solution is a global minimizer for h^ < h^ < has o,nd the normal 
solution is a global minimizer for /iq > has ■ 

In S3 n A2 the symmetric solution is locally stable for < ho < h and 
locally unstable for h < /ig < hs, the asymmetric solution is locally stable 
fof has < ho < has o-^d the normal solution is locally stable for ho > has ■ 
Moreover, there exists he £ {hs,h) such that the upper symmetric solution 
is a global minimizer for < ho < he, the asymmetric solution is a global 
minimizer for he < ho < has cLnd the normal solution is a global minimizer 
for ho > has- 

This conjecture is consistent with physical observations of superheating 
and supercoohng. The critical field h is often called superheating field be- 
cause when starting from a symmetric superconducting solution near h = 0, 
the material will remain in this locally stable state until h is reached. Then 
as ho is increased above h, the maximum of / has a discontinuity and bulk 
superconductivity is destroyed. 

Now, on the contrary, starting from normal solution, when ho is de- 
creased from infinity, the material remains in the normal state until the 
nucleation field {hg or has) is reached. If the branch of bifurcating solutions 
is unstable, this field is called supercooling field because nucleation occurs 
through a jump in the maximum of /: the material gets into the symmet- 
ric state which is the global minimizer. These two phenomena give rise to 
hysteresis loops. 



For all the cases they have studied, Bolley-Helffer [13| have established 
the local stability or instability of the branches bifurcating from normal 
solutions. In fact, when they proved the bifucation to be subcritical, they 
also proved local instability of the branch, and stability for supercritical 
branches, except in ^3 n A2 and Si Ci A2. Note that in ^3 n A2 and Si n 
A2, though the bifurcation is supercritical, they proved that the symmetric 
solution is unstable. 

For a fixed and ho < ^o(^) defined earlier in Section 3.5.2, Aftalion 
[Q] proved that the global minimizer tends to a symmetric solution when 
K tends to infinity. Indeed, in this range corresponding to Si D A2, there 
are no asymmetric solutions for ho < hQ{a), but for /iq > hQ{a), it is an 
interesting open problem to prove that both symmetric and asymmetric 
solutions coexist and the asymmetric solution is the global minimizer. 
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As we pointed out in Section 4, in 5*3 n A2 and 82^1 A2, the formation 
of sheaths and vortices is related to the existence of stable asymmetric so- 
lutions. The Hastings- Troy result ^^ gives a first step in the mathematical 
analysis of this phenomenon: Hastings and Troy ^^ gave a first step in the 
mathematical analysis of this phenomenon since they proved that for fixed 
K > I/V2.OI and large a, there is a range of h values for which the asym- 
metric solution is the global minimizer. As we have said before in Section 
4, when the asymmetric solution is a global minimizer, it is not known for 
which regime of (a, k), it gives rise to surface sheath or to vortices. 

We have provided a comprehensive study of the structure of symmetric 
and asymmetric solutions for the entire range < a, k < 00. Throughout 
we have tied together our new observations with previously known results. 
At the core of our new results are the existence of the triple point and the 
quintuple point. These points provide the birth of families of asymmetric 
and symmetric solutions via bifurcation. Already, one of the authors, to- 
gether with S.J. Chapman, has begun further investigations of the structure 
of both of these points and this has resulted in two new papers y], ||3[. 

Acknowledgement. The authors are very grateful to P.J.McKenna for 
mentionning the existence of AUTO. 
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